Synchronized clusters in coupled map networks: Stability analysis 
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We study self-organized (s-) and driven (d-) synchronization in coupled map networks for some 
simple networks, namely two and three node networks and their natural generalization to globally 
coupled and complete bipartite networks. We use both linear stability analysis and Lyapunov 
function approach for this study and determine stability conditions for synchronization. We see 
that most of the features of coupled dynamics of small networks with two or three nodes, are carried 
over to the larger networks of the same type. The phase diagrams for the networks studied here 
have features very similar to the different kinds of networks studied in Ref. [T^|. The analysis of 
the dynamics of the difference variable corresponding to any two nodes shows that when the two 
nodes are in driven synchronization, all the coupling terms cancel out whereas when they are in 
self-organized synchronization, the direct coupling term between the two nodes adds an extra decay 
term while the other couplings cancel out. 
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I. INTRODUCTION 

Recently it has been observed that several complex sys- 
tems have underlying structures that are described by 
networks or graphs having interesting properties Q, 0- 
Many of these naturally occurring large and complex net- 
works come under some universal classes and they can be 
simulated with simple mathematical models, viz small- 
world networks Q, scale- free networks Q etc. These 
models are based on simple physical considerations and 
have attracted a lot of attention from physics community 
as they give simple algorithms to generate graphs which 
resemble actual networks found in many diverse systems 
such as the nervous systems social groups |(| .world 
wide web 0, metabolic networks 8], food webs and 
citation networks [Toj . 

Several networks in real world consist of dynamical el- 
ements interacting with each other and they have large 
number of degrees of freedom. Synchronization in dy- 
namical systems with many degrees of freedom has at- 
tracted much attention in recent decades |l5l 0. In 
Refs. 0, 0] two of us have presented detailed analy- 
sis and numerical results of phase synchronization and 
cluster formation in coupled maps on different networks. 
Starting from random initial conditions the asymptotic 
behaviour of these coupled map networks (CMNs) has re- 
vealed that there are two different mechanisms leading to 
two types of synchronized clusters. First there are clus- 
ters with dominant intra-cluster couplings which are re- 
ferred as self- organized (s-) synchronization and secondly 
there are clusters with dominant inter-cluster coupling 
which are referred as driven (d-) synchronization |lllll2fl. 
The numerical studies reveal several clusters with both 



types as well as clusters of mixed type where both mech- 
anisms contribute. There are also situations where ideal 
clusters of both types are observed. An analysis of sim- 
ple networks with two and three nodes indicates that the 
self-organized behaviour has its origin in the decay term 
arising due to intra-cluster couplings in the dynamics of 
the difference variables while the driven behaviour has its 
origin in the cancellation of the inter-cluster couplings in 
the dynamics of the difference variables. 

In the present paper we study the dynamics of some 
simple networks analytically and numerically with a view 
to get a better understanding of the two mechanisms of 
cluster formation discussed above. Mainly we study the 
asymptotic stability of s- and d-synchronization in net- 
works with small number of nodes, i.e. two and three 
nodes, and extension of these small networks to networks 
with large number of nodes. As an example of large net- 
works showings-synchronized clusters we take globally 
coupled maps |13| and for d-synchronized clusters we take 
complete bipartite coupled maps |14| . 

The paper is organized as follows. In Section II, we 
state the model and discuss general features of the sta- 
bility of synchronized states. Section III discusses two 
networks showing s-synchronization, i.e. two nodes net- 
work and globally connected network. Section IV consid- 
ers two networks showing d-synchronization, i.e. three 
nodes bipartite network and complete bipartite network. 
Section V concludes the paper. 



II. MODEL OF A COUPLED MAP NETWORK 
(CMN) 
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Consider a network of N nodes and N c connections (or 
couplings) between the nodes. Let each node of the net- 
work be assigned a dynamical variable x 1 , i = 1,2, ... ,N. 
The dynamical evolution of coupled maps can be written 



2 



as m 



L t+i 



= (l-e)f(xt) + 



N 
3 = 1 



(1) 



where x\ is the dynamical variable of the i-th node at the 
i-th time step, e is the coupling strength (0 < e < 1), C 
is the adjacency matrix with elements Cij taking values 
1 or depending upon whether i and j are connected 
or not. C is a symmetric matrix with diagonal elements 
zero, hi = ^2 Cij is the degree of node i. The function 
f(x) defines the local nonlinear map and the function 
g{x) defines the nature of coupling between the nodes. 
For most of our results we use the logistic map 

f(x) = fj,x(l - x) 

for illustration. We use two types of the coupling func- 
tions, 



X 



(2a) 
(2b) 



We refer to the first type of coupling as linear and the 
latter as nonlinear. 

Synchronization: Synchronization of two dynamical 
variables or systems is indicated by the appearance of 
some relation between functionals of two processes due 
to interaction IH 13 El El) In this paper we will mostly 
concentrate on exact synchronization, where the values of 
the dynamical variables associated with nodes are equal. 
We can see that the fully synchronized state of a network, 
x\ = x\ = ■ ■■ = x N = Xt, is a solution to Eq. (Q. This 
fully synchronized state lies along a one-dimensional di- 
agonal in the N dimensional phase space of the dynamical 
variables. We define a synchronized cluster as a cluster 
of nodes in which all pairs of nodes are synchronized. As 
stated in the introduction, we can identify two different 
mechanism of cluster formation in CMNs. Self-organized 
synchronization which leads to clusters with dominant 
intra-cluster couplings and d-synchronization which leads 
to clusters with dominant inter-cluster couplings 11]. 

In this paper we study synchronization in some simple 
networks. We use two types of analysis to determine 
the stability of synchronized state. First is the linear 
stability analysis [17lll8lll9ll20l.l2ll . l22l] and the second is 
Lyapunov functional |22l|2j,|2J. We now briefly discuss 
these two methods. 

Linear Stability Analysis: The evolution of tangent 



vector Sz t = [5x1, 6a 
be written as, 



, 5x 



NlT 
t J ' 



along a trajectory can 



Sz t 



+i 



JtSzt 



where Jt is the Jacobian matrix at time t, (Jt)ij — 
dx\ +l l 'dx\. If J t is a diagonal matrix or the similarity 
transformation which diagonalizes Jt is independent of t, 



then the Lyapunov exponents can be written in terms of 
the eigenvalues of Jt as follows, 



A; 



lim -VlnlA^)! 

r— »oo t — * 



(3) 



t=i 



where Ai(t) is the i-th eigenvalue of Jacobian matrix at 
time t. If Jt does not satisfy the above mentioned condi- 
tions then it is necessary to consider product of Jacobian 
matrices to obtain Lyapunov exponents. 

To study the stability of synchronized state it is 
sufficient to consider transverse Lyapunov exponents 
which characterize the behavior of infinitesimal vectors 
transversal to the synchronized manifold, and these de- 
termine the stability of a synchronized state [2^. If all 
the transverse Lyapunov exponents are negative then the 
synchronized state is stable. 

Global Stability Analysis: Condition for global stability 
can be derived using the Lyapunov functional methods 
[231. |24| ] . Global stability in a neighborhood of an equilib- 
rium (stable) point is confirmed if there exist a positive 
definite function defined in that neighborhood, whose to- 
tal time derivative is negative semi-definite. 

To get the conditions for the global stability of syn- 
chronization of two trajectories x\ and x 3 t , we define Lya- 
punov function as, 

(4) 



JV2 



Clearly Vij(t) > and the equality holds only when the 
nodes i and j are exactly synchronized. For the asymp- 
totic global stability of the synchronized state, Lyapunov 
function should satisfy the following condition in the re- 
gion of stability, 

V t+1 < V t 

This condition can also be written as, 
V t+1 



< 1. 



(5) 



III. SELF-ORGANIZED SYNCHRONIZATION 

We first consider the simplest and smallest network 
showing s-cluster, i.e synchronization of two coupled 
nodes which is obviously of self-organized type. As a 
generalization of two nodes network to large networks, 
we study the s-synchronization in globally coupled maps. 

A. Coupled network with N = 2 

We begin by taking the simplest case where number 
of nodes is two and these two nodes are coupled with 
each other E3 ■ The dynamics of the two nodes can be 
rewritten as (Eq. Q), 



b t+i 



b t+i 



(1 
(1 



e)f(xl) 
e)f(x 2 t ) 



tg(x t ) 
tg(x\) 



(G) 
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1. Linear stability analysis 

Following Ref. 0, we first define addition and dif- 
ference variables as follows, 



Arf < 0. Thus the range of stability of the synchronized 
state is given by 



1 



e~ A " 1 + e 
— < e < 0.5 — 



(11) 



St 



d t = 



(7) 



Dynamical evolution for these newly defined variables is 
given by 



1 



st+i = 



H+l 



[f(st + dt) + f(st-ck)} 
dt) +g(st - d t )] 



-[f( St + d t )-f(s t -dt) 
[g(s t + d t ) -g(s t -d t )] 



(8a) 



For synchronous orbits to be observed, the fully synchro- 
nized state d t — i.e. x\ = x\ = x t — s t , should be a 
stable attractor. The Jacobian matrix for the synchro- 
nized state is, 



Jt 



(l-e)f(x t ) + eg'(x t ) 






(1 - e)f(x t ) - eg'(x t ) 



where the prime indicates the derivative of the function. 
The above Jacobian is a diagonal matrix and Lyapunov 
exponents can be easily written in terms of eigenvalues of 
product of such Jacobian matrices calculated at different 
time. The two Lyapunov exponents are 

X s , d = Urn -Y / ln\(l-e)f(xt)±eg'(xt)\ (9) 

t=l 

The synchronous orbits are stable if Lyapunov exponent 
corresponding to the difference variable d tl i.e. Xd or the 
transverse Lyapunov exponent, is negative. If the other 
Lyapunov exponent A s is positive then the synchronous 
orbits are chaotic while if it is negative then they are 
periodic. 

Coupling function g(x) = f(x): Two coupled maps with 
g(x) — f(x) type of coupling are studied extensively in 
the literature both analytically and numerically 0, [25| . 
For the synchronized state Lyapunov exponent A s is 
nothing but the Lyapunov exponent for uncoupled lo- 
gistic map (A u ) and the other Lyapunov exponent can 
be written in terms of the A„ and from Eq. @ we get 



A., 



A,/ 



A,, 



lim 

r — >oo T 



-X>l/'(* t )l 



t=l 

ln|l-2e|+A„ 



(10a) 
(10b) 



The synchronous orbits are stable if Lyapunov exponent 
corresponding to the difference variable dt is negative, i.e. 



For logistic map with (j, = 4, this gives 0.25 < e < 0.75 
as the range for the stability of the synchronized state. 
Coupling function g(x) — x: For g(x) — x type of 
coupling, numerical results show that as the coupled 
nodes evolve, dynamics shows different types of synchro- 
nized and periodic behaviors depending upon the cou- 
pling strength e and the parameter of the map /. First 
let us start with the general case where coupled dynamics 
lies on a synchronized attractor. Using Eq. @, Lyapunov 
exponents can be easily written as 



lim -y>|(l-e)/'(s t )±e| 



r — >oo 7" * 

For stable synchronous orbits Lyapunov exponent corre- 
sponding to the difference variable, i.e. A^, should be 
negative. Now we consider some special cases when cou- 
pled dynamics lies on periodic or fixed point attractor. 

Here, as well as for other networks considered in this 
paper, we will restrict ourselves to period two orbits. 
Higher periodic orbits exist but are difficult to treat an- 
alytically. Also major features of phase diagram are un- 
derstood by using fixed point and period two orbits. 

Case I. Synchronization to period two orbit: Consider 
the special case where the solution of Eq. © or Eq. lj%a)l 
is a periodic orbit of period two, with the difference vari- 
able given by dt+\ = d t = and addition variable given 
by s t+ 2 — s t = si, s t+ 3 = s t+ i = s 2 . Eigenvalues for 
the product matrix J\ J2, where J\ and J2 are Jacobian 
matrices for consecutive time steps, are given by 



Ai, 



(1 - t?Af-2 + £ 2 ±(1 - e)e(f[ + }' 2 ) (12) 



where f[ and f 2 are derivatives of / at the two peri- 
odic points (si,g? = 0) and (s 2 ,d = 0) respectively. The 
e range for which the dynamical evolution gives a sta- 
ble periodic orbit, is obtained when the modulus of the 
eigenvalues for matrix J\J 2 are less than one. For local 
dynamics given by logistic map f(x) — fj,x(l — x), the 
two periodic points are, 



S1.2 



For 



1 + e + n(l - e) ± v /e(e-2)( A i-l) 2 -3-2/z + ^ 2 



2 M (1 - e) 



Sl,2 = 



5 - 3e + V5 - 18e 



9e 2 



8(1 - e) 



(13) 



(14) 



which gives the coupling strength range 0.18.. < e < 
0.24.. for which the periodic orbit, (s 1; d = 0) and (s 2 , d — 
0), is stable. 
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Case II. Period two orbit: There is a range of e values 
that give the following stable period two behaviour, 



1 2 _ 1 

x t — X t+1 — x t+2 
T 2 _ 1 2 _ vl 



b t+2 



Lyapunov exponents for this periodic state can be found 
from the eigenvalues of the product of Jacobians at the 
two periodic points. Jacobian matrix at (Xf, X%) is given 
by, 



(15) 



where /{ and f 2 are the derivative of / at X\ and X\ 
respectively. Eigenvalues of the product matrix Ji J 2 are, 



A ll2 =(c±(l-e)v70t) S 



(16) 



If f[f 2 < which is the interesting case, then the condi- 
tion for the stability of the periodic orbit become 



M - 1 



< e < 1. 



(17) 



For f(x) — (jlx{1 — a;), Xf and X\ which satisfy Eq. i(T5|) 
are given by 



A l,2 



(l + 7T)±V( 1 + ^ 2 -|( 1 + ^ 



(18) 



For stable periodic orbits modulus of both eigenvalues 
Ai and A2 are less than one which gives the coupling 
strength range for stability as, 



If the expression in the square bracket on the RHS is less 
than one then the synchronized state is stable. 

For the nonlinear coupling function g(x) = f(x) the 
expression (I21[l for Lyapunov function simplifies and we 
get 



X^ x~ 



= (1 - 2 £ ) 2 [/'(,^/"(^) + 0[{x\ - x 2 ) 2 ] 2 

(22) 

If the expression in the square bracket on the RHS is 
bounded then there will always some range of e values 
around 0.5 for which the synchronized state will be sta- 
ble. For logistic map fix) = (ix(l — x), we get, 

V t+1 = V t (l-2e) 2 (i 2 [\-{x\+x 2 )] 2 

Using < x\ + x\ < 2, we get the following range of e 
values for which the synchronization condition given by 
Eq. JSJ is satisfied, 



(23) 



For (j, = 4, it gives the coupling strength range 0.325 < 
e < 0.675. However, a better e range can be obtained by 
putting more realistic bounds for X^ ~\- Xj- cLS 

where X = sup t (|l — x\ — x 2 \). 



B. Globally coupled networks 



1 



[i 2 - 2 t i - 3 



< e < 1 



(19) 



For n — 4, we get coupling strength range (0.6 < e < 1) 
for which periodic orbit as given in Eq. i|15|) is stable and 
coupled dynamics lies on a periodic attractor. 



2. Lyapunov function analysis 

From Eqs. (@J and ©, Lyapunov function for two 
nodes is written as 



Globally coupled networks have all pairs of nodes con- 
nected to each other i.e. N c = N(N — l)/2 where N is 
the number of nodes. For such global coupling we write 
our dynamical model as 



xl +1 = {l-e)f{x\) + w — i ]T 9(4) (25) 



Let the state x\ = x 2 = ■ ■ ■ = x^ = x t be the fully 
synchronized state. We now consider the stability of this 
state. 



V t % = v t+1 



= [(1 - e)(/04) - fix*)) - e(g(x$) - g(x 2 )^0) 

Using Taylor expansion of f{x\) and g(x\ ) about x 2 , we 
get Lyapunov function at time t + 1 as, 

V t+1 = V t [(l-e)f'(x 2 )-eg'(x 2 ) 



+ ^r^((l-e)f"(x 2 )-e)g"(x 2 )) 



+ 0[{x)-x 2 r 



(21) 



1. Linear Stability Analysis 

Jacobian matrix at time t for the fully synchronized 
state is 



Ji 



({l-e)f t jfrgi ... jfrgl \ 
*M (l-e)/ t ' ... T^g't 



V w=l9t w=i9t ■■■ (l-e)/ t '/ 
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where f' t and g' t are the derivative at the synchronous 
value Xf. Eigenvectors of the above Jacobian matrix are, 



E m = [ exp[2m—), exp(^m — ) exp{2NTTi — ) 



where m = 0, 1, . . . , N — 1 and T denotes the transpose. 
From these eigenvectors we find that Jt has an eigenvalue 
(1 — e)fl + eg' t and (N — l)-fold degenerate eigenvalues 
(1 — ejfj. — N e _ 1 g'f Lyapunov exponents can be written 
in terms of eigenvalues of Jacobian matrix as, 



Ai 
A 2 



lim i]rin|(l-e)/; + 6 9 ; 

r — >oc T L — * 

t=l 

A3 • • • = Xn 
lim — In 

r — >oo T — ' 



N - 1 



9t 



(26a) 



(26b) 



Lyapunov exponents A2, A3 . . . A n are transverse Lya- 
punov exponents since they characterize the behaviour 
of infinitesimal vectors transverse to the synchroniza- 
tion manifold. For stability of the synchronous orbits 
all transverse Lyapunov exponents should be negative. 

Coupling Function g{x) = f{x): Globally coupled 
maps are studied extensively with g(x) — f(x) type of 
coupling and linear stability analysis is done to decide 
the stability of fully synchronized solution [TtI HH , l2fj| . 
Lyapunov exponents are given by, 



Ai = A u = lim y2hx\f'(x t )\ 



At — X-> 



A at = X u + ln(l 



N 



N - 1 



(27a) 



e) (27b) 



Except X u other Lyapunov exponents are transverse to 
the synchronization manifold. The critical value of cou- 
pling strength e c beyond which the synchronized state 
with all nodes synchronized with each other is stable, is 
given by, 



N-l 

N 



(28) 



For large N with f(x) = 4x(l — x), we have X u = In 2. 
So from above expression we get e c = 0.5. 

Coupling function g(x) — x: From the expressions 
(I26all and (|26bp . it is difficult to determine when the 
synchronous orbits are stable. Rather it is easier to de- 
termine the stability of the synchronous orbits using Lya- 
punov function which we will consider in the next subsec- 
tion. Here, we consider some special cases when coupled 
dynamics of the fully synchronized state lies on periodic 
or fixed point attractors. 

Case I. Synchronization to fixed point: First we con- 
sider the fixed point X* — x t of the fully synchronized 
state. The fixed point is given by X* — f(X*). The con- 
ditions for the stability of the synchronous fixed point 



are In |(1 - e)f + e\ < and In |(1 - ejf - ^| < 0. 
where /' is the derivative of / at X*. For /' < which 
is the interesting case, the fixed point is easily found to 
be stable in the range 



1/1-1 
1 

N-l 



l/'l 



< £ < 1. 



(29) 



For N = 2 this solution is not stable and the range of sta- 
bility increases with N. This is a surprising result since 
with increasing N the number of transverse eigenvectors 
along which the synchronized solution can become un- 
stable also increases. This feature of increasing range of 
stability with N is more general and will be noticed for 
other solutions also as will be discussed in subsequent 
analysis. For large N synchronous fixed point is stable 
for 



in 1 
i/'i 



< e < 1 



(30) 



For f(x) given by logistic map, X* = 1 — 1/fi, f'(X*) = 
2 — [i and range of stability of the synchronous fixed point 
is given by 



[i — 3 



M 2 N 1 _ 1 



< e < 1. 



(31) 



Case II. Synchronization to period two orbit: We now 
consider synchronous period two solution. The Lyapunov 
exponents can be obtained using Eqs. I)26al) and l|26b() . 
The conditions for the stability of the period two solution 
are 



lrx|((l-e)/{ + e )((l-e)^ + e)| <<32h) 



In 



N ■ 



1) O 1 -^ 



N - 1 



<32b) 



For f{x) given by logistic map, the period two synchro- 
nized solution is the same as given by Eqs. Ijl3|l and 114(1 . 
For \i = 4, the range of stability of this solution is 



1 - 



10 



< e < 



1 

N-l 



60 



30 
N-l 



{N-l) 2 



2 

N-l 



(33) 



(N-l) 1 



For N — 2, this gives the e range (0.18. .,0.24..) for the 
stability of period two solution as noted in Casel of two 
node case with g(x) = x. For large N, this range of 
stability for e values expands to (0.18. .,0.28..) pt| . 



C. Lyapunov Function Analysis 

From Eqs. 10} and (|25p. we write Lyapunov function 
for any two nodes of globally coupled network as, 



v t+i 



(f(xl)-f(a 
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Performing Taylor expansion around x\ , we get 
Vt+i 



V 



(±-z)f(4) 



N- 1 



g'(4) 



+ ■ 



2 

0[(4 



(l-e)f"{x{)- 



N - 1 



g"(4) 



„2\21 



(34) 



Coupling function g(x) = f(x): In this case the ex- 
pression l|34|l simplifies to 



Vi 



t+i 



v t 



= 1- 



N 



N - 1 



f'(4) 



-f"(4) 



0[{x\ - x 2 t f 



If the expression in the square bracket on the RHS is 
bounded then for large N there will be a critical value 
of e beyond which the condition JSJ will be satisfied and 
the globally synchronized state will be stable. For f(x) = 
fj,x(l—x) and using (0 < x\+x\ < 2), we get the following 
range of coupling strength values for which the globally 
synchronized state is stable. 



N-l 

N 



1 



< e < 1< 



N - 1 



N 



1 



(35) 



For fi — 4 and for very large N, coupling strength range 
is 0.75 < e < 1. A better e range is obtained by putting 



more realistic bounds as X < x\ 
gives the range of stability as, 



< X + , which 



IV. DRIVEN SYNCHRONIZATION 

Driven synchronization leads to clusters with dominant 
inter-cluster couplings. For the ideal d-clusters, there are 
only inter-cluster connections with no connections be- 
tween the constituents of the same cluster ^lj- A com- 
plete bipartite network consists of two sets of nodes with 
each node of one set connected with all the nodes of the 
other set. Clearly this type of network is an ideal exam- 
ple for studying d-synchronization. We take a bipartite 
network consisting of two sets of nodes, N — Nl + N2, 
with each node of Nl connected to every node of N2 and 
there are no connections between the nodes of the same 
set ^ij. We study dynamics of coupled maps on such 
type of bipartite network and determine the stability cri- 
teria for formation of d-synchronized clusters. Our model 
for coupled complete bipartite network can be written as, 



H+l 



N 



^-e)f(x\) + — J2 9(4), 

j=Nl+l 

fori = I,-- - ,N1 

Nl 

{l-e)f{x\) + ^£9(4), 

3=1 

fori = JV1 + 1,- •• ,N. 



(39) 



where all terms are having the same meaning as defined 
for Eq. 0. 



Network with three nodes 



N-l 

N 



1 



1 



< e < 1 



(36) 



where A = max(|l - X+\, |1 - X~\). 

Coupling function g(x) — x: In this case, for logistic 
map, the expression l|34|) simplifies to 



v t+i 
V t 



(l-e)(l-(xl+xl))- 



N- 1 



Since < x\ -I- x\ < 2 we get a range of e values for which 
the globally synchronized state is stable (Eq. |J5}) 

fi — 1 ^ _ ^ M + 1 



i 

N-l 



< e < 



l 

N-l 



which for large N reduces to 



< e < 



(37) 



(38) 



For fj, = 4, we get the coupling strength range 0.75 < 
e < 1 for which the globally synchronized state is sta- 
ble. We also note that for TV = 2 the condition l|37|l for 
synchronization is not satisfied for any value of e < 1. 



First we take a simple network of three nodes with 
Nl = 2, N2 = 1 which is the smallest possible network to 
show the behaviour displayed by Eq. I|39|l . The evolution 
equations can be written as 



x t+i 
x t+i 

3 



X 



t+1 



(1 
(1 
(1 



*)/(*?) 
e)f(x 3 t ) 



eg (4) 

f (9(4) 



9(4)) 



(40) 



Here, the d-synchronized state corresponds to the nodes 
1 and 2 being synchronized with each other but not with 
node 3. 



1. Linear Stability Analysis 

Here, we use addition and difference variables St and 
dt as defined in Eq. (0 for the first two nodes. Thus 
Eqs. Ij40|l can be rewritten as, 



st+i = 



1-e 
2 

1 -e 



[/(*t + dt)+/(«t-dt)] + eff(s?)(41a) 



d t +i = —[/(«t + dt)-/(«t-dt)] (41b) 
4+i = (l-t)f(4) + ^[9(st + d t )+g(s t -d t ^lc) 
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Jacobian matrix for the d-synchronized state (st,dt 
0,Xf) is given by, 

(1 - eg' 2 

■h = | (1 - e)f[ j (42) 
eg[ (1 - e)/ 2 

f[ and g[ are derivatives at s t = x t = x\ = x\ and 
/2 and g' 2 are derivatives at {x\\ Lyapunov exponent 
corresponding to the difference variable is found as, 



where \ r is, 



X d = ln(l-e) + A c 



= lim j>|/'(« t )| 



(43) 



(44) 



The coupling strength range for which the d-synchronized 
solution is stable, l.e A^ < 0, is given by, 



1 



f 



exp A f 



< e 



(45) 



Coupling function g{x) = f(x): We have investigated 
the network containing three nodes numerically with 
g{x) — f(x). For logistic map with /i = 4, we find that 
nodes get d-synchronized for coupling strength ranges 
0.15 < e < 0.2 and 0.5 < e < 1. When we investi- 
gate further these two e ranges we find that in the lower 
range the behaviour is mostly periodic while in the upper 
range it is periodic in the middle portion and is chaotic 
at both ends. For 0.5 < e < 1 it is easy to see why the 
synchronized dynamics gives a stable attractor. Since 
ln(l - e) < - In 2 for e > 0.5, from Eq. we see that 
till A c is less than In 2, the value of Lyapunov exponent 
for an isolated logistic map at /i = 4, the synchronized so- 
lution will be stable. The stability of the d-synchronized 



state for 0.15 < e < 0.2 appears to be because of either 
periodic attractor or chaotic attractor close the periodic 
attractor with very small value of A c . 

Now we discuss the special cases where nodes are syn- 
chronized with variables showing periodic or fixed point 
behaviour. For this study we use the original variables. 

Case I. Synchronization to Fixed point: Nodes get syn- 
chronized to a fixed point with one set of nodes having 
one value and the other set of nodes having a different 
value such as (x\ = x\ — X*) and (x^ — X 2 ). Eigenval- 
ues of Jacobian matrix iB^ with g(x) = f(x) are, 



Ai 

A 2 , 3 



^(/i+/a) 
1 



(46a) 



4(1 - 2e)/(/£46b) 



where /{ and f 2 are the derivatives at X\ and X 2 respec- 
tively. By putting the condition that the magnitude of 
the above eigenvalues should be less than one, we get the 
e range for which the fixed point state is stable. When 
the first two nodes synchronize with each other, the three 
coupled maps system behaves like just two coupled maps 
and within this e range all the solutions are those of the 
two coupled maps. The expressions for X\ and X 2 are 



_ (1 - /i + 2^e) ± y/(l - /i + 2 M e) 2 - 4e(l - /i + 2//e) | 

(47) 

From Eqs. (|4(ja|l and l|46b|) the coupling strength values 
for which the fixed point solution is stable must satisfy 



|(l-e)/{|<l 



and the following condition 




< e < 



3 + 2n(n - 2) + t/(3 + 2n(n - 2)) 2 - Afijfi - 2)(fi - if 
MP - 2) 



(48) 



r 



Case II. Synchronization to period two orbit: Consider 
a periodic orbit of period two where the first two nodes 
are d-synchronized and have the value Xf , and the third 
node has has the value X 2 . The period two orbit is ob- 
tained by interchanging these two values at successive 
time steps. The Jacobian matrix for this periodic orbit 
can be written as J = J\ J 2 where 



(1 - e)f[ ef 2 
Ji = I (1 - e)f{ eft 



(49) 



where /{ and f 2 are derivatives at X\ and X 2 respec- 
tively and J 2 is obtained by interchanging suffixes 1 and 
2 in the expression for J\ . Eigenvalues of J are given by 

Ai = (l-e) 2 f[f 2 (50a) 

\ 2 



A 2 , 3 = ^(/i + /2)±Vj(/(+^) 2 + (l-2 e )»|) 
For f{x) given by logistic map, the periodic points are 
X\ 2 = a ± (51) 




FIG. 1: Phase space diagram showing different features of coupled dynamics in the two parameter space of 11 and e for three 
nodes bipartite network with logistic map as local map and coupling function g(x) = f[x). Different regions are T. Turbulent 
region, DP. Driven periodic, DF. Driven fixed point, DC. Driven Chaotic, S. Self organized region and F. Fixed point. Region 
boundaries are determined based on the asymptotic behaviour using several initial conditions, synchronization behaviour and 
the largest Lyapunov exponent. The dashed lines indicate uncertainties in determining the boundaries. The inset shows the 
phase diagram for the entire range of parameter fi from to 4. 



where 



1 / 1 

S = - 1 + — ; : 

2 V M(l " 2e) 



After imposing the condition l-A.1,2,3 1 < 1 we can get the 
e range for which the period two orbit is stable. The first 
eigenvalue gives us the condition that 



1- e < 



II — 2e] 



+n- 2e^)(4e - 1 + \i - 2efi) ' 



(52) 



The other two eigenvalues give us the same condition as 
in Eq. I|48[) except that e in this inequality is substituted 
by 1 — e similar to the case of two coupled maps [25j ■ For 
logistic map with fj, = 4, this range is given by 0.16.. < 
e<0.20... 

Case III. Synchronization of all three nodes: It is also 
possible that all the three nodes get synchronized (s- 
synchronization) . The eigenvalues of Jacobian matrix 
for this synchronized state (x\ = x\ = x\ = Xt) can 
simply be written from Eqs. (|46a[) and (|46 bl . by putting 
f{ = fi = f'{ x t)i which gives Lyapunov exponents as, 

1 T 

Ai = X u = - lim V] In I /'(a t ) I (53a) 

t=l 

A 2 = ln(l-e)+A„ (53b) 
A 3 = ln(l - 2e) + A„ (53c) 

where X u is the Lyapunov exponent for an uncoupled 
map f(x). For synchronous orbits to be stable, the two 



Lyapunov exponents (A2 and A3) corresponding to the 
transverse eigenvectors should be negative. We get the 
following e range for which both the transverse eigenval- 
ues are negative, 



<£< 2 (1 



(54) 



In this region coupled dynamics may lie on a chaotic or 
a periodic attractor depending upon the value of A u . For 
f(x) = fix(l—x) with fj, = 4, A lt = In 2 and thus the range 
of stability of the s-synchronization of all three nodes is 
(0.5, 0.75) and the dynamics is chaotic. 

Phase diagram in ji — e plane: Fig. ^ shows differ- 
ent phases in the fJ- — e plane for three nodes bipartite 
network with g{x) — f(x). For /1 < 3 we get a fixed 
point solution. To understand the remaining phase di- 
agram consider the line ijl = 4. Fig. [21 shows two sets 
of differences between the values of variables, Ix 1 — x 2 \ 



(open circles) and ja; 1 



x 



(crosses) as a function of the 



coupling strength e. Bipartite d-synchronized state and 
global s-synchronized state are clearly seen. Fig. OJa) 
shows largest Lyapunov exponent and Fig. 0fa) shows 
the fractions of inter- and intra- couplings, j 'inter and 
f intra, a function of e for mu = 4. The two quantities 
/intra and /inter act as measures for the intra-cluster cou- 
plings and the inter-cluster couplings and are defined as 
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FIG. 2: The figure shows the variation of two sets of differ- 



(open circles) and 



xf\ (crosses) for a 



three nodes network as a function of the coupling strength e 
for f(x) = nx(l — x) with /i = 4 and g(x) = f(x). For each e, 
100 values of the differences are plotted after an initial tran- 
sient. 



0111, 



/intra 
fi liter 



N c 

•Winter 



N c 



(55a) 
(55b) 



where -ZVj n t ra and -ZVj ntcr are the numbers of intra- and 
inter-cluster couplings respectively. In iVi n ter ; couplings 
between two isolated nodes are not included. Initially 
for small coupling strength values nodes are in turbulent 
region with no cluster formation at all (region T). As the 
coupling strength increases beyond a critical e c we get 
bipartite d-synchronized state (region DP). This region 
corresponds to the Case II of period two orbits discussed 
above in this subsection. When the coupling strength 
increases further we get a reappearance of turbulent re- 
gion. As the coupling strength is increased further all the 
nodes are synchronized giving one cluster with global s- 
synchronization (region S). This region corresponds to 
the case III discussed above and the range of e for this 
region is given by Eq. I|54|) . In this region the coupled 
dynamics lies on a chaotic attractor. In the last three 



regions (DC, DP and DC), we get driven bipartite syn- 
chronization. The middle region (DP) displays driven 
fixed point solution (case I) while the other regions show 
chaotic behaviour. 

For [i < 4, the coupling strength range where all 
the nodes are synchronized (region S) increases in size 
and the coupling strength ranges where nodes show d- 
synchronization shift towards the nearest boundaries e = 
and e = 1. 

Coupling function g(x) — x: Numerical analysis shows 
periodic solutions which we now consider in detail. When 
the first two nodes are d-synchronized the solutions of the 
dynamical equations are similar to the two nodes case. 
This simplifies the analysis considerably. 

Case I. Driven synchronization to period two orbit: 
First we consider the case where the coupled dynamics 
shows periodic attractor with period two behaviour and 
the variables take the values 



(56) 



Using the product of Jacobians of Eq. (|42|l for two con- 
secutive time steps the eigenvalues for this periodic orbit 
can be easily obtained. The eigenvalue A^ associated 
with the difference variable d t is given by 



A d = (1 - eff[f 2 , 



(57) 



The other two eigenvalues are the eigenvalues of product 
matrix J1J2, where J\ is given by Eq. l|15|l and hence 
the eigenvalues are given by Eq. (|16|) . The solution of 
the periodic orbit is the same as for two coupled maps 
(Eq. (|18(l ). Using all the three eigenvalues of Jacobian 
matrix, the coupling strength range, for which the peri- 
odic orbit given by Eq. (|56[) is stable, is given by 



max 1 



1 



\m 



1 



1 



< e < 1 



(58) 



The periodic points X\ and X2 are also the periodic 
points of uncoupled map / (Eq. 1)13(1 ). For logistic map, 
we get the following range of coupling strength in terms 
of logistic map parameter /j,, 



1 



< 6 < 1 



(59) 



which is the same as Eq. (|19|) . The lower bound of e 
matches exactly with the boundary between the regions 
IV-DP and IV-DQ of Fig. 1 of Ref. [13 and regions DP 
and DQ of Fig. 

Case II. Synchronization of all three nodes: All three 
nodes get synchronized for a small coupling strength re- 
gion with dynamics lying on periodic orbits of period two 
such that 



u t+i 



u t+i 



4 



X\ 



(60) 



10 



0.8 



0.6 



0.4 



0.2 



0.2 



1 1 




(a) 


0.8 


— 


1 1 

1 
1 


(b) - 


. T i 




d : 


0.6 


1 — 

• 
• 


, j 

TTT — n 

Ill 1 J • TT T 1 — 1 


• IV — L) 


\ ! ° ; t / 


S 






• 
• 


iv-p • IV s 


X* 


'•.! ! y 




* 






1 

\! 


• • 

* 

• 






• 
• 


0.4 




t ;% 


• 
• 






• 


A 






• 


- 1 i 






0.2 




i ' 





























j - -. 




.•! 




• 

• 






. * • . - 

• 


• 

• 








-0.2 








1 1 










i.i I 





0.2 0.4 0.6 0.8 1 



0.2 



0.4 0.6 

C 



0.8 



FIG. 3: (a) The figure shows the largest Lyapunov exponent A as a function of the coupling strength e for the three nodes 
bipartite network with logistic map as local map with fi = 4 and g(x) = f(x). (b) Same as for (a) but for a bipartite network 
of 50 nodes with Nl = N2 = 25. 



The eigenvalue of the Jacobian for this periodic orbit for 
the difference variable is simply (1 — t) 2 f\f 2 and the other 
two eigenvalues are the same as for two coupled maps and 
are given by Eq. (fT^)) . The periodic points are given by 
Eq. (|14|) . The coupling strength range for which all three 
nodes are synchronized for logistic map with fi = 4, is 
0.18... < e < 0.24... which is same as Case I for two 
coupled maps with g{x) = x. 

Phase diagram in n — e space: Fig. shows differ- 
ent phases in the fi — e plane for three nodes bipartite 
network with g(x) = x. For /i < 3 we get a fixed 
point solution. To understand the remaining phase di- 
agram consider the line /j, — 4. Fig. shows two sets 
of differences between the values of variables, {x 1 — x 2 \ 
(open circles) and {x 1 — x 3 \ (crosses) as a function of the 
coupling strength e. Bipartite d-synchronized state and 
global s-synchronized state are clearly seen. Fig. Efa) 
shows largest Lyapunov exponent and Fig. ffia) shows 
the fractions of inter- and intra- couplings, f inter and 
f intra, as a function of e for mu = 4. Initially for small 
coupling strength values nodes are in turbulent region 
with no cluster formation at all (region T). As the cou- 
pling strength increases beyond a critical e c , we get global 
s-state (region S). This region is the case II considered 



above. When the coupling strength increases further we 
get a reappearance of turbulent region. In the last two 
regions we get driven bipartite synchronization (regions 
DQ and DP). The last region corresponds to case I of 
period two discussed above in this subsection and the 
critical coupling strength for it is given by Eq. ||SSJ. 

For /i < 4, the coupling strength region for d- 
synchronization gets wider and that for s-synchronization 
gets thiner with a shift towards e = 0. 



B. Complete bipartite networks 

Let us now consider a complete bipartite network of 
N = Nl + N2 nodes and dynamics defined by Eq. 
We define a bipartite synchronized state of the bipartite 
network as the one that has that all Nl elements of the 
first set synchronized to some value, say X\(i), and all 
N2 elements of the second set synchronized to some other 
value, say X2(t). Linear stability analysis of the bipar- 
tite synchronized state can be done using the Jacobian 
matrix, 
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FIG. 4: (a) The figure shows the fractions of inter- and intra-cluster couplings, f inter and f intra, as a function of the coupling 
strength e for the three nodes bipartite network, logistic map as local map with fi = 4 and g(x) — f(x). The values are obtained 
by averaging over 50 random initial conditions, (b) Same as for (a) but for a bipartite network of 50 nodes with Nl — N2 — 25. 
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(61) 



■■■ (l-e)/^ 



where and g 2 are the derivative of g(x) at Xl and X 2 eigenvalues of the above Jacobian matrix can be divided 
respectively. It is easy to see that the eigenvectors and into three sets, A, B, C, as, 

I 



set Eigenvectors Eigenvalue No of eigenvalues condition 

A (ai . . . ajvi,0- • -0) (1 - e)f[ Nl - 1 Sa = 

B {0...0,{h...M £/3 = 
C (a,...a,0,...0) - 2 



r 



Here, (ai, . . . ajvi), (/?i, • • -/3jv2) and a, /3 are complex column. The two eigenvalues corresponding to the set C 
numbers satisfying the conditions specified in the last 
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FIG. 5: Phase space diagram showing different features of 
coupled dynamics in the two parameter space of /x and e for 
three nodes bipartite network with logistic map as local map 
and coupling function g(x) = x. Different regions are T. Tur- 
bulent region, DP. Driven periodic, DF. Driven fixed point, 
DQ. Driven Quasiperiodic, DC. Driven Chaotic, S. Self or- 
ganized region and F. Fixed point. Region boundaries are 
determined based on the asymptotic behaviour using several 
initial conditions, synchronization behaviour and the largest 
Lyapunov exponent. The dashed lines indicate uncertainties 
in determining the boundaries. The inset shows the phase 
diagram for the entire range of parameter fi from to 4. 



are the eigenvalues of the matrix 

(1 - 0/f ^ 
effi (1 " e)& 



(62) 



We make the following observations, (a) The three sets 
of eigenvectors, A, B, C, are orthogonal to each other and 
the total space of eigenvectors may be written as a direct 
sum of these three sets, (b) The three sets of eigenvec- 
tors do not mix with each other under time evolution, 
(c) The eigenvectors belonging to the first two sets (A 
and B) are also the eigenvectors of the product of any 
number of Jacobian matrices under time evolution, (d) 
The two sets of eigenvectors A and B are transverse to 
the synchronization manifold which is defined by the last 
set of eigenvectors C. 

Lyapunov exponents corresponding to the transverse 
eigenvectors can be easily written as 



A; 



Xi = 



ln|(l-e)| + i lim £ln|/'(*i)l, 



t=l 

fori = 1, ... , Nl — 1 



ln|(l-e)| + i lim f>|/'(X 2 )|, 



fori = Nl, 



t=i 
,N-2. 



(63a) 



(63b) 
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FIG. 6: The figure shows the variation of two sets of dif- 



ference, 



(open circles) and \x t — xf\ (crosses) for a 



three nodes network as a function of the coupling strength e 
for f(x) = nx(l — x) with fi — 4 and g(x) = x. For each e, 100 
values of the differences are plotted after an initial transient. 



The synchronized state is stable provided the transverse 
Lyapunov exponents are negative. If /' is bounded as 
is the case for the logistic map then from Eqs. (|63|) we 
see that for e larger than some critical value, e c (< 1), 
bipartite synchronized state will be stable. Note that this 
bipartite synchronized state will be stable even if one or 
both the remaining Lyapunov exponents corresponding 
to the set C are positive, i.e. the trajectories are chaotic. 

Now we study the periodic behaviour of coupled dy- 
namics of Eq. (|39|l . Here we consider two cases, fixed 
point attractor and periodic attractor with period two, 
as we studied for three nodes. 



1. Coupling function g(x) — f(x) 

The fixed point bipartite synchronized solution with 
one set of nodes taking value XI and the other set X| , is 
given by Eq. i|47|) ■ Jacobian matrix for this solution gives 
three sets of eigenvalues, first set of Nl — 1 degenerate 
eigenvalues (1 — e)f{, second set of N2 — 1 degenerate 
eigenvalues (1 — ejf^ and third set of two eigenvalues 
given by Eq. (|46b(l . The conditions for the stability of 
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FIG. 7: (a) The figure shows the largest Lyapunov exponent A as a function of the coupling strength e for the three nodes 
bipartite network with logistic map as local map with /j, — 4 and g(x) = x. (b) Same as for (a) but for a bipartite network of 
50 nodes with Nl = N2 = 25. 



this solution is given by Eq. (gSJ) an d 1(1 - 2 | < 1. 

The bipartite synchronous period two solution is ob- 
tained when one set of nodes take the value Xf and the 
other set of nodes take the value X% and the two values 
alternate in time. This solution is the same as given by 
Eq. Tne eigenvalues (Eqs. (|50ajl and (|50b|) . with 

the N — 2 fold degeneracy of Ai) and the conditions for 
stability of the solution are the same as for Case II of 
three nodes with g(x) = f(x). 

The global synchronized solution where all nodes of the 
bipartite network are synchronized, has the same Lya- 
punov exponents (except the degeneracy) as in Eqs. (|53cl) 
and the stability criterion is again given by Eq. I|54|). 



2. Coupling function g(x) = x 

Consider a periodic orbit of period two with bipartite 
synchronized state where all nodes of one set take value 
X\ and and the other set X% and two values alternate 
in time. The solutions is same as for case I of three 
nodes bipartite network. The stability ranges are given 
by Eqs. {55)1 and Pty. 



C. Lyapunov Function Analysis 

For the complete bipartite networks, Lyapunov func- 
tion as defined by Eq. Q, for any two nodes belonging 
to the same set is given by 



v t+l — 



[{l- ( ){f(x l t )-f{xl)t 



(64) 



Expanding around x\ gives the ratio of Lyapunov func- 
tions at two successive times, 



V 13 
v t+i 



f'(4) 



l ±f"{4) + o{{x\-4f 



If the term in the square bracket on the RHS is bounded 
then there will be a critical value of e beyond which 

yij 

-^tj- < 1 and thus the bipartite synchronized state will 
be stable. 

We see that for d-synchronization, Lyapunov function 
for any pair of nodes does not depend on the size of the 
complete bipartite network and type of the coupling be- 
cause in the expression for Lyapunov function, contribu- 
tion of such couplings cancel out. This is not the case for 
globally coupled networks where contribution of the cou- 
pling terms for the two nodes under consideration do not 
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FIG. 8: (a) The figure shows the fractions of inter- and intra-cluster couplings, f inter and /intra, as a function of the coupling 
strength e for the three nodes bipartite network, logistic map as local map with /i = 4 and g(x) = x. The values are obtained 
by averaging over 50 random initial conditions, (b) Same as for (a) but for a bipartite network of 50 nodes with Nl = N2 — 25. 



exactly cancel and the size of the network has an effect 
on the asymptotic behaviour. 

For the logistic map and using < x\ + x\ < 2 and 
Eq. we get the following range for e values for syn- 
chronization of nodes i and j, 

^—1 < e < 1 (65) 
M 

A better e range can be obtained by taking a more ap- 
propriate boundary as X~ < x\ + x\ < X + , which gives 

1 

1 < e < 1 

fiA 

where A is defined after Eq. Q36|l . 

D. Journey from small network to large network 

To understand the features of d-synchronization in bi- 
partite network and comparison of this large network 
with the network of three nodes, here we briefly present 
the numerical results for logistic map with fi = 4. 
Fig. I^Jb) plots the largest Lyapunov exponent as function 



of e and Fig. EJb) plots f inter and f intra as a function of 
e for a bipartite network with g(x) — f(x) and jj, = 4. 
Comparing with the three nodes network (Figs. G3a) 
and^Ja)) we find that the major difference is in the range 
0.20.. < e < 0.5 where the turbulent region reappears for 
the three nodes network while the mixed region having 
both d- and s-clusters is observed for the larger bipar- 
tite networks. The region boundaries for the period two, 
globally s-state and fixed point are the same for the three 
nodes and larger bipartite networks and are the same as 
given for cases II, III and I respectively for three nodes 
bipartite network with g(x) = f(x). We see that most of 
the dynamical features and synchronized cluster forma- 
tion in coupled map bipartite networks with three nodes 
get carried over to the larger bipartite networks. 

As in the case g(x) = f(x) there is a similarity between 
the three nodes bipartite network and the larger bipar- 
tite networks for g(x) = x. Fig. E{b) plots the largest 
Lyapunov exponent as a function of e and Fig.|SJb) plots 
finter and f intra as a function of e for a bipartite net- 
work with g(x) = x and /i = 4. Comparing with the 
three nodes network (Figs. Efa) and Ufa)) we find that 
the main difference is in the range 0.24.. < e < 0.35.. 
where the turbulent region reappears for the three nodes 
network while d-synchronization appears for the larger 
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bipartite networks. The other regions show a very sim- 
ilar behaviour for both three nodes and larger bipartite 
networks. 

E. S- and D-synchronization 

The analysis presented so far shades some light on 
the dynamical origin of the two types of synchronization 
namely self-organized and driven, that we have studied. 
In s-synchronization the clusters have dominant intra- 
cluster couplings while in d-synchronization they have 
dominant inter-cluster couplings. We consider two nodes 
and globally coupled networks as simple examples of s- 
synchronization while three nodes and complete bipartite 
networks as simple examples of d-synchronization. The 
dynamics of the difference variable and Lyapunov func- 
tion analysis are useful for understanding the difference 
between the two mechanisms of cluster formation. 

From Eq. l|8a|) (b) we see that in the dynamics of the 
difference variable dt for the two nodes network the cou- 
pling term adds an extra decay term. This is also seen 
from the expression l(21|l for Lyapunov function. On the 
other hand, from Eq. I|41b(l we see that in the dynamics 
of the difference variable for the three nodes network the 
coupling terms with the third variable cancel out (see also 
Eq. 164(1 for Lyapunov function). The situation is more 
complicated when we consider larger networks. The d- 
synchronization shows the same trend i.e. cancellation 
of the coupling terms in the dynamics of the difference 
variables and as well as in the expression for Lyapunov 
function (Eq. JMJ). On the other hand, Eq. lj54*j) for Lya- 
punov function for the globally coupled maps shows that 
the direct coupling term between the two nodes under 
consideration adds an extra term in the difference vari- 
able while the coupling terms to other variables cancel 
out. 

The analysis presented in this paper is for exact syn- 
chronization while Ref. ^3 considers phase synchronized 



clusters. However, we feel that the dynamical origin for 
the two mechanisms of cluster formation should be sim- 
ilar in both cases. This is supported by the very similar 
features of plots of phase space, largest Lyapunov expo- 
nent, f inter and f intra for the three nodes and complete 
bipartite networks and the corresponding plots for sev- 
eral networks considered in Ref. 1121. 



V. CONCLUSION 

We have studied s- and d-synchronization in coupled 
map networks using some simple networks, namely two 
and three node networks and their natural generaliza- 
tion to globally coupled and bipartite networks respec- 
tively. For this study we use both linear stability anal- 
ysis and Lyapunov function approach. We see that 
for the difference variable for any two nodes that are 
in d-synchronization all the coupling terms cancel out 
whereas when they are in s-synchronization though cou- 
pling terms for couplings to other nodes may cancel, the 
coupling terms corresponding to the direct coupling be- 
tween the two nodes under consideration do not cancel. 

The phase diagrams for three nodes network has fea- 
tures very similar to the different kinds of networks stud- 
ied in Ref. [T^ . (The phase diagrams for two nodes, not 
shown in this paper, have similar features except the d- 
synchronization which needs at least three nodes to man- 
ifest itself.) Under the time evolution, the coupled dy- 
namics lies on different periodic or chaotic attractors with 
varying coupling strength values. The type of coupling 
function plays an important role in the time evolution 
of the coupled dynamics. We also see that most of the 
features of coupled dynamics of small networks with two 
or three nodes, are carried over to the larger networks of 
the same type. 
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